Mechanical state, material properties and continuous description of an epithelial tissue 
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During development, epithelial tissues undergo extensive morphogenesis based on coordinated 
changes of cell shape and position over time. Continuum mechanics describes tissue mechanical 
state and shape changes in terms of strain and stress. It accounts for individual cell properties 
using only a few spatially averaged material parameters. To determine the mechanical state and 
parameters in the Drosophila pupa dorsal thorax epithelium, we sever in vivo the adherens junctions 
around a disk-shaped domain comprising typically hundred cells. This enables a direct measurement 
of the strain along different orientations at once. The amplitude and anisotropy of the strain increase 
during development. We also measure the stress to viscosity ratio and similarly find an increase 
in amplitude and anisotropy. The relaxation time is of order of ten seconds. We propose a space- 
time, continuous model of the relaxation. Good agreement with experimental data validates the 
description of the epithelial domain as a continuous, linear, visco-elastic material. We discuss 
the relevant time and length scales. Another material parameter, the ratio of external friction to 
internal viscosity, is estimated by fitting the initial velocity profile. Together, our results contribute 
to quantify forces and displacements, and their time evolution during morphogenesis. 
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I. INTRODUCTION 

An epithelial tissue is a sheet of cells that acts as a bar- 
rier, separating for instance the outside and the inside of 
a multicellular organism. Its biological function relies in 
part on the formation of a network of adherens junction 
belts, connected to the acto-myosin cytoskeleton, where 
cells adhere to each other [1], and which transmits me- 
chanical information over several cell diameters [1-3]. A 
key issue is to understand and model the role of tissue 
mechanics (forces, displacements, and their time evolu- 
tion) in the coordinated changes of cell shape and posi- 
tion which determine morphogenetic flows at the tissue 
level [1, 4, 5]. 

Several models describe tissues using continuum me- 
chanics [3, 6-12]. One precondition is the existence of a 
mesoscopic scale defining a domain over which averages 
of cell properties are well-defined [6, 13]. This description 
further relies on the assumption that the tissue mechani- 
cal state can be quantified, at the same mesoscopic scale, 
by two variables [6]: the stress characterizing in which 
directions, and to what extent, the domain is under ten- 
sion or under compression; and the strain characterizing 
how far the domain geometry is from that of a relaxed 
state. Such description has the advantage of accounting 
for individual cell properties using only a few spatially 
averaged parameters, which determine for instance how 
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fast and through which succession of states the domain 
reacts to external solicitations. 

Beyond simple observation, in vivo mechanical mea- 
surement techniques include elastography [14], photoc- 
lasticity [15], magnetic micromanipulation [16], tonome- 
try [17] or nanoindentation [18]. A large literature (for 
review, see e.g. [19]) has established laser ablation of in- 
dividual cell junctions as a tool to measure the tensions 
within an epithelium, in particular during Drosophila 
dorsal closure [3, 20-24]. This technique has allowed to 
measure the material relaxation time t and contributed 
to a better characterization of morphogenetic processes 
in Drosophila [25-28]. 

The analysis of single cell junction ablation is usually 
based on reasonable assumptions [19, 20]: that the tissue 
is at mechanical equilibrium just before severing; that the 
ablation is effective in removing at least part of the cell 
junction tension; and that during relaxation the velocity 
remains small (Reynolds number much smaller than 1). 
Within these assumptions, the initial retraction velocity 
yields the value of the cell junction tension removed by 
the ablation, up to an unknown prefactor which depends 
on the dissipation [25-29]. 

The tensor a denotes the component of the stress re- 
moved by ablation. It arises as an average over several 
individual cell junctions within a region of the tissue. By 
ablating straight lines in either of two perpendicular di- 
rections [30, 31] or by performing statistics on single cell 
junction ablations in several directions [3], it is possi- 
ble to measure step by step the anisotropy of the ratio 
a/rj (where n denotes the tissue viscosity). The stress 
to viscosity ratio a/n is the initial strain rate after the 
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FIG. 1. Large-size annular severing in a fly dorsal thorax. 

(a) Developmental stages of a fruit fly Drosophila 
melanogaster. Top: larva. Middle: pupa, with pupal 
case removed (ESM Figure 1) and cuticle kept intact. Bot- 
tom: adult. Dashed lines represent the midline (symmetry 
axis). The X axis is antero-posterior: anterior (head) 
towards the left, posterior (abdomen) towards the right. 
The Y axis is medio-lateral. Crosses: approximate positions 
of severing, along the midline in the scutellum (blue) and 
off-axis in the scutum (black). 

(b) Epithelial cell apical junctions marked by E- 
CadherimGFP just before the severing, t = 0~ , here in 
an old pupa (see text for classification). Blue circles: two 
concentric circles define the annular severed region; the 
distance between circles corresponds to about one cell size. 
Yellow arrows: macrochaete used as spatial references to 
position the severed region. 

(c) First image after severing, t = 1 s. Yellow: fitted ellipse 
[36] (see Methods), (d) Time t = 30 s after the severing, 
showing a larger opening along Y than X. Bars = 1 mm (a), 
10 /im (b-d). 



severing: the rate of displacement that the epithelium 
would spontaneously undergo if it were free. 

The continuous description applies to other cellular 
materials such as foams [32, 33]. In the present work, 
we investigate experimentally whether it also applies to 
epithelial tissues. More specifically, we ask whether the 
epithelium strain, stress and material properties can be 
directly measured on a mesoscopic scale. 

We address these questions in a model system, the dor- 
sal thorax epithelium of Drosophila pupa. The pupa is 
the life stage during which the larva starves, feeds upon 
its reserves and metamorphoses into an adult [34]. The 
dorsal thorax is a single layer of cells whose apical sur- 
faces, surrounded by the adherens junction belts, face the 
cuticle, which protects the pupa; their basal surfaces face 



the hemolymph, which acts as a nutrient transporter. It 
is composed of a large anterior region, the scutum, and 
a posterior tip, the scutellum. The thorax, and espe- 
cially the scutellum, undergoes extensive morphogenetic 
changes during the metamorphosis that shapes the adult 
fly (Figure 1 (a)), making it a useful model to decipher 
the mechanisms that control tissue morphogenesis dur- 
ing development [29]. It has a bilateral symmetry axis: 
the midline (axis X, dashed line in Figure 1 (a-b)). It 
displays large cells, the macrochaete (precursor cells of 
adult sensory hairs, Figure 1 (b)), whose positions are 
precisely reproducible [35] . 



II. ANNULAR SEVERING EXPERIMENTS 

We introduce an original type of severing experiment. 
We sever by short pulse laser the adherens junctions in 
an annular region around a ^30 /jm radius circular tissue 
domain (see Figure 1 (b-d), Methods and Movies 1-3). 
The inner and outer wound margins display a comparable 
speed, displacement and anisotropy. Note that the same 
anisotropy in speed and displacement results in the outer 
ellipse boundary having minor and major axes orthogonal 
to the minor and major axes of the inner one. After 
severing, the wound heals within tens of minutes and the 
pupa develops to adult. 

Experiments are done in the scutellum (blue cross in 
Figure 1 (a)). We classify them into three groups, ac- 
cording to the developmental age at which the severing is 
performed: earlier than 15-16 h, approximately 18-20 h, 
and later than 24-26 h after pupa formation. For brevity 
we call them "young" (N = 8), "middle-aged" (N = 5) 
and "old" (N = 10) pupae, respectively. Similar exper- 
iments can be performed in other locations. As a proof 
of principle, we include three experiments performed in 
the scutum (black cross in Figure 1 (a)) of middle-aged 
(N = 2) and old (N = 1) pupae. 

We analyse the retraction of the inner tissue domain 
and its margin, within the set of assumptions relevant 
for the retraction of a single severed cell-junction. The 
boundary of the retracting domain is fitted at each time 
point by an ellipse [36] (see Methods and Figure 1 (c-d)), 
yielding values of the X and Y axes as a function of time 
(Figure 2). 

Our protocol offers several advantages, (i) The ex- 
periment directly yields measurements averaged over the 
severed tissue domain, comprising typically a hundred of 
cells. As discussed below, the domain size is thus an ad- 
equate mesoscopic scale, crossing-over from detailed cell- 
level to tissue-level continuous descriptions [6, 13]. (ii) 
The comparison between the initial and final states yields 
a direct measurement of the strain which existed before 
severing, (iii) Several cell-cell junctions in various direc- 
tions are severed at once (rather than step by step) and 
in the same pupa. Hence one experiment yields a direct 
measurement of all stress to viscosity tensor components. 
It takes only a few tens of seconds, the typical value of 
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FIG. 2. Model-independent measurement of strain e, stress to 
viscosity ratio cr/r/, and relaxation time r. The wound margin 
position (ellipse semi-axis) is plotted versus time after sever- 
ing; data from Figure 1 (b-d), ESM Movie 1, along the Y axis. 
The difference between the initial (L) and final (Ly) positions 
(blue dashed lines) directly yields the value of e = ln(L/Ly ). 
The initial velocity — ^jr(t = 0) (slope of the orange dashed 
line) yields an estimate of La/rj. Inset: Velocity, estimated 
by finite differences of successive positions, versus the posi- 
tion during the first 30 s. An arrow indicates the direction of 
increasing time t. The slope of a linear fit (purple line) yields 
the inverse of the relaxation time, 1/r (see Equation 2). 



the relaxation time, (iv) Last but not least, by isolating 
in time and space an epithelium domain from its neigh- 
boring cells, we obtain a complete set of well-defined ini- 
tial, boundary and final conditions for a spatio-temporal 
model of the retracting domain. We thus introduce a 
continuous model, in the spirit of Mayer et al. [30] and 
taking into account a coupling between space and time 
dependences. As explained below, fitting such a model 
to the displacements in the domain bulk tests the contin- 
uous description, suggests an interpretation of the relax- 
ation time, and probes the competition between internal 
viscosity and external friction. 



III. RESULTS 

A. Model-independent measurement of strain, 
stress to viscosity ratio and relaxation time 

Severing the epithelium reveals the displacement that 
the epithelium would spontaneously undergo if it were 
free. We measure the strain using a logarithmic definition 
[32, 37, 38], Hencky's "true strain" [39]: 
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where L is the initial radius, and is the ellipse semi- 
axis at the end of the relaxation. Similarly £yy = 



ln(L/L^?)- The advantages of this definition are that 
[32, 33]: it is valid at all amplitudes; it is adapted for ten- 
sorial measurements; and finally it increases the range of 
validity of the linear elasticity approximation which ap- 
plies here even to the highest strain value we measure, 
0.45. 

Each severing is followed by a relaxation to a final do- 
main strictly smaller than the initial disk (Figure 1 (d)), 
indicating that before the severing the tissue had a pos- 
itive strain in all directions. Since the midline is a sym- 
metry axis, we expect that the strain axes are parallel 
and perpendicular to it. We check that this is the case 
for the fitted ellipse axes, and that accordingly the shear 
strain sxy is indistinguishable from 0. We thus plot exx 
and Syy (Figure 3 (a)), obtained with absolute precision 
better than 10 -2 . The values are clustered according 
to the three pupa age groups: e is low and isotropic at 
young age (green cluster), moderate and isotropic at mid- 
dle age (red), high and anisotropic (eyy > Exx) at old 
age (blue). 

Similarly, in each experiment, the initial retraction of 
the wound margins (Figure 1 (c)) indicates the sign of the 
stress: before the severing the tissue was under tensile 
stress in all directions. This is reminiscent of positive 
cell junction tensions observed in the wing [25] and in the 
notum [29]. The initial retraction velocity divided by L 
yields the stress to viscosity ratio a/rj (Figure 2), in a way 
which is likely to be independent of any tissue rheological 
model, as suggested by the following argument. When 
ablating a single cell-cell junction, force balance shows 
that the force F exerted by the ablated junction on the 
neighbouring vertex is proportional to the initial recoil 
velocity v. F = — jv, where 7 is a friction coefficient 
[19]. Coarse-graining this relation, and assuming that the 
main source of dissipation is identical in our experiments 
which sever many junctions at once, yields a ~ r]v/L, 
where the factor L is introduced by dimensional analysis. 
Results for the stress to viscosity ratio a/rj (Figure 3 (b)) 
are similar to those for e: its amplitude and anisotropy 
are initially small and increase with age; the values are 
clustered according to the three age groups; no shear 
stress is detected. Values of a/77 span almost two decades, 
reflecting the sensitivity and precision of the method. 

The relaxation time, t, is obtained from the time evo- 
lution of the ellipse axes. Instead of fitting an exponential 
to the data, we plot the velocity versus position (Inset of 
Figure 2) . This method is robust and independent of any 
assumption. A linear regression of the velocity as a func- 
tion of position yields a slope of —1/r, as seen, e.g., in 
the equation: 
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We find that the relaxation time is approximately 
isotropic (ESM Figure 3 (A)) and use (t x + t y )/2 as 
an estimate of r. It is close to 10 s and only slightly 
varies with the pupa age group (Figure 3 (c), vertical 
axis) . 
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FIG. 3. Mechanical state and material properties. Color 
code according to pupa development ages: green, young; red, 
middle-aged; blue, old. The edges of the rectangular regions 
represent the mean values ± standard deviations for each 
group. The experiments in the scutum are in yellow, (a) 
Strain anisotropy: eyy versus Exx- Note the difference in 
horizontal and vertical scales; the solid line is the first bisec- 
trix Y — X, indicating the reference for isotropy. (b) Same 
for the severed stress to viscosity ratio ojr\. (c) Relaxation 
time r and dimensionless friction to viscosity ratio £; values 
are the averages of the measurements along the X and Y axes 
(ESM Figure 3). Note the semi-log scale. The blue rectangle 
takes into account two very small values of £, of order 10 -3 
and 10 -4 (below the plotted range). 



B. Space-time model and measurement of the 
friction to viscosity ratio 

To describe the bulk of the retracting domain, we track 
the displacements of around a hundred features in a band 
of tissue (Figure 4 and Methods). Spatially averaging 
their positions by binning them into 8 equal-size bins 
improves the signal to noise ratio. We obtain the initial 
velocity profile versus position (Figure 5). For old pupae, 
it is spatially linear; for middle-aged pupae it is spatially 
nonlinear; in the young pupae, immediately after severing 
only the boundaries move significantly. 

To account for these observations, we formulate a 
spatio-temporal, viscoelastic, Kelvin- Voigt model in one 
dimension of space, where the coordinate z denotes either 
X or Y, and t denotes time after severing (see ESM). 
Thanks to the definition (1), the strain e(z,t) is related 
without approximation to the velocity v(z,t) through: 
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We model the effect of friction of the epithelium against 
the hemolymph and the cuticle as an external fluid fric- 
tion with coefficient £. We find that when internal viscos- 
ity dominates, the strain remains uniform and all parts of 
the tissue relax exponentially with the same viscoelastic 
relaxation time r = r\j ' E where E is the Young modulus 
(ESM Figure 4 (A)). When external friction dominates, 
strain diffuses from the boundaries with a diffusion co- 
efficient D = -EyC and remains thus inhomogeneous at 
times t smaller than t d = L 2 /D = (L 2 /E (ESM Fig- 
ure 4 (B)). In the general case, the dynamical equation 
for the strain field s(z,t) reads: 
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with the conditions e(z,t = 0) = Sq initially, e(z = 
±L,t) = £oexp(— i/r) at boundaries, and e(z,t — > oo) = 
at the end. The relative importance of external friction 
versus internal viscosity is quantified by the dimension- 
less parameter: 
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Solving Equation 4 at short time and integrating over 
space Equation 3 at t = yields the initial velocity pro- 
file: 
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whose curvature is controlled by £. We thus estimate £ by 
fitting a hyperbolic sine to experimental initial velocity 
profiles (Figure 5). As for the relaxation time, we check 
that the values £x and £y determined in the directions 
X and Y are roughly consistent (ESM Figure 3 (B)): we 
thus use (£x + £y)/2 as an estimate of £. 
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FIG. 4. Tracking of feature positions, (a) Positions of features (blue squares), old pupa, data of Figure 1 (c) after image 
denoising. The logarithms of intensity levels are represented in grey scale, and contrast is inverted for clarity, (b) Tracking for 
the first 30 s (time color-coded from blue to red), inside a rectangle along the Y axis, (c) Same for the X axis. 



We find that £ varies more than r (Figure 3 (c), hori- 
zontal axis). The middle aged group is characterized by 
intermediate values of £ (Figure 5, red), of order of 10. 
According to Equation 5, this indicates a possible effect 
of the external friction on the relaxation for scales larger 
than ~ 10 fxm. For both other age groups, the 

values of £ should be considered as bounds on the order 
of magnitude. For old pupae we find £ smaller than a few 
times unity (Figure 5, blue). For young pupae, we find £ 
larger than a few tens or hundred (Figure 5, green). 

This model depends on three parameters measured in- 
dependently of each other: the initial strain eq, the vis- 
coelastic time r, and the friction to viscosity ratio £. To 
validate the model, we use it to numerically compute the 



whole space-time map of strains and displacements. We 
simulate each experiment using the values of r and £ 
measured as described above. Figure 6 shows that we 
find a good agreement with experimental data. More- 
over, we plot the square deviation between experimental 
and numerical data as a function of the parameter values 
(r, £) (see Methods). We check that the error landscapes 
we obtain are broadly consistent with our experimental 
estimates of £ and r (ESM Figure 5). 
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FIG. 5. Model-dependent measurement of the friction to vis- 
cosity ratio £. Initial velocity profiles immediately after sev- 
ering are plotted versus initial position prior to severing, for 
three typical pupae: young (green), middle-aged (red) and 
old (blue). Open squares: features, from Figure 4 (b). Closed 
triangles: same, spatially averaged in 8 bins. Lines: fit by a 
sinh function (Equation 6), yielding for £ a value either above, 
in, or below the measurable range: £ > 60 (green), £ ~ 14±5 
(red), £ < 0.5 (blue), respectively. 



agrees with values provided by single cell junction abla- 
tion: 10-100 s for the relaxation times, 0.1 to 1 /mi.s -1 
for the initial velocities in other Drosophila tissues [3, 25- 
28]. The model allows to interpret r as a viscoelastic 
time. Its order of magnitude is typical of the viscoelastic 
times corresponding to cell internal degrees of freedom 
[40], also measured in cell aggregates [41, 42]. 

Conversely external degrees of freedom (such as cell- 
cell rearrangements or cell divisions) would have much 
larger viscoelastic times: in cell aggregates they are typ- 
ically of order of several hours [11, 41, 43]. The time 
scales of morphogenesis in this tissue are also of order of 
hours [29], as observed in other tissues as well, such as 
the Drosophila pupal wing [12]. Indeed continuous mod- 
els of morphogenetic processes presuppose the existence 
of such a separation of time scales. 

The present experimental data are well described as- 
suming that a single viscoelastic time is relevant. Other 
experiments display a wide distribution of viscoelastic 
times [3, 28]: the model can easily be modified so as to 
take this feature into account. Note that in the pres- 
ence of friction the model already predicts that different 
parts of the severed domain relax at different rates (ESM 
Figure 4 (B)). 
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Measurements of e, cr/77, r and £ are performed as av- 
erages over the severed disk. The requirements on the 
averaging scale are that it should be (i) large enough to 
obtain a sufficient signal to noise ratio yielding significant 
results, and (ii) small enough to be able to focus simulta- 
neously on the whole severed ring despite the epithelium 
surface curvature which is observed in old pupae. The 
Y — > — Y symmetry with respect to the midline observed 
in Figures 1 (b-d), 5, and 6 acts as a control of the sig- 
nal to noise ratio. The experimentally available range 
of scales corresponds to 60-100 cell areas, in other words 
the semi-axis length L corresponds to 4-5 cell diameters. 
For old pupae, scales smaller than L are also probed, cor- 
responding to the scale of bins (at least L/2 or L/A, see 
Figs. 5 and 6). 



FIG. 6. Validation of the model. Blue triangles: positions of 
features from Figure 4 (b), spatially averaged in bins versus 
time after severing; bars: standard deviation. Red lines: nu- 
merical resolution of Equation 3, using the measured values 
£ = 0.08, r = 15.8 s. 



IV. DISCUSSION 

A. Time scale 

We report relaxation times r in the range 5 to 20 s, 
and initial velocities in the range 0.1 to 0.6 /^m.s -1 . This 



C. Friction 

We experimentally determine the order of magnitude 
of the dimensionless friction to viscosity ratio, £ (Figure 3 
(c), horizontal axis). Obtaining more precise values of £ 
would depend on model ingredients, such as depth of ab- 
lation along the apico-basal axis, detailed cell geometry, 
tissue compressibility, or type of friction. 

For simplicity we have modeled by a fluid friction 
[30] the dissipation on either surfaces of the epithelium, 
namely the hemolymph at basal side and the cuticle at 
apical side. We do not observe any threshold effect due 
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to solid friction, and the agreement between model and 
data justifies a posteriori our assumption. 

For young pupae, the initial velocity field is not signif- 
icantly higher than the noise level (Figure 5, green). In 
this age group, it cannot be used to validate the model, 
and the estimate of £ strongly depends on the motion of 
domain parts located close to the boundary (Equation 6). 

D. Continuous description 

In the dorsal thorax, the severed domain fulfills the 
theoretical requirements for a continuous description: it 
contains a number of cells much larger than one (large 
enough to allow for averaging), while being much smaller 
than the total number of cells in the whole tissue. 

One could ask whether averages at scale L are mean- 
ingful, due to variation of cell properties from cell to cell. 
We speculate that if heterogeneities within the domain 
dominated the average, the domain boundary could in 
principle adopt any arbitrary shape compatible with the 
Y — > —Y symmetry. However, in practice, we observe 
that the domain boundary at the beginning, as well as at 
the end, of the relaxation also displays a X — > —X sym- 
metry (Figure 1 (b-d) and Movies 1-3) and can be fitted 
by an ellipse: this is consistent with the measurements 
of a/ rj and e as tensors averaged over scale L. 

E. Model 

Limits of our model lie in its simplifications: one 
dimensional treatment, linear elasticity and viscosity, 
isotropy and homogeneity of material parameters includ- 
ing friction, identical properties of all cells contained in 
the severed region. However, the agreement between 
model and experiment is already good, and validates 
these approximations a posteriori. We emphasize that 
one major goal of this model is to provide an order of 
magnitude of the friction to viscosity ratio £. We do not 
claim that this model is unique, nor that it is general. In 
the future we hope to further test the model ingredients 
and relevance on this epithelium and other quasi-planar 
ones. 

Since its ingredients are general, we expect our model 
to be of broad relevance when studying the relaxation 
of quasi-planar epithelia over similar time and length 
scales. However, the simplifications involved deserve fur- 
ther comments. In fact, completely different hypotheses 
could also be invoked to explain the spatial non-linearity 
in the velocity profile observed for young and middle-aged 
pupae. For instance, the viscous component of biologi- 
cal tissues could be a non-linear function of time, strain, 
and strain rate [6]. A strain-dependent viscosity coeffi- 
cient may lead to a non-linear velocity profile without 
the need to invoke any external friction. Alternatively, 
a spatially heterogeneous friction might also explain the 
observed velocity profiles. Since a linear velocity profile 



is observed in older pupae, the amplitude of these addi- 
tional, more complex ingredients would have to decrease 
strongly during development. 

We believe that our model has the merit of simplic- 
ity. It successfully accounts for the description of the 
velocity profiles at any time. Each of its hypotheses 
could be tested experimentally and, if experimental re- 
sults require it, hypotheses could be relaxed to lead to 
more detailed explorations. We expect that additional 
ingredients could strongly affect the value of £, but that 
other measurements presented here (namely e, a/rj and 
r) should be robust. The model thus provides a flexible 
and robust route to a mechanical description of quasi- 
planar epithelia. 



F. Perspectives 

Building upon classical single-junction laser ablations, 
our experiments enable to measure the time evolution 
of the mechanical state (strain, and stress to viscosity 
ratio) and material properties (relaxation time, and fric- 
tion to viscosity ratio) during pupa metamorphosis, as 
well as their anisotropy. This provides relevant ingredi- 
ents for the modeling of a tissue as a continuous, linear, 
viscoelastic material. 

Moreover, our approach provides a method to ex- 
plore the biological causality of the observed properties, 
by analysing mutant phenotypes or by micro-injecting 
drugs. We conjecture that r would be affected by low 
doses of drugs known to modify cytoskeletal rheology 
such as Cytochalasin D [44], CK-666 [45],Y-27632 [46] 
which respectively inhibit actin polymerization, Arp2/3 
complex activity or Rho-Kinase activity. We anticipate 
that modulating cell-cell adhesion by knock-down of E- 
Cadherin or Catenins [47] could affect the viscosity. Al- 
though we cannot modulate the cell-cuticle adhesion, we 
could tune cell-matrix adhesion by knock-down of inte- 
grin or its linker to the cytoskelcton [47] to affect the 
friction and thereby modify £. 

In principle, by establishing a complete map of a at 
different positions, one could determine the length scale 
at which a varies within the tissue, and check a posteriori 
whether it is larger than the scale L used to measure a. 
However, the present experiments measure the stress up 
to a dissipative prefactor, the tissue viscosity. Similarly, 
single cell-junction ablation experiments measure force 
up to a prefactor (a friction coefficient), rarely discussed 
or measured in the literature. To separate the stress and 
viscosity variations, alternative methods to directly mea- 
sure forces and viscosities in live epithelia are thus called 
for. These include the measurement of cell junction ten- 
sions from movies observations [48], where the unknown 
prefactor, namely the average tension over the whole im- 
age, is by definition the same for all junctions. External 
mechanical manipulation [16-18] yields direct determi- 
nation of out-of-plane elastic and/or viscous moduli, but 
the in-plane moduli are only indirectly determined. An- 
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other possible strategy could rely on microrhcology [49- 
51]: either passive, by tracking the brownian diffusion of 
beads within the cells; or active, by moving (magneti- 
cally or optically) a bead back and forth within a cell at 
a specified frequency. Such measurements could yield a 
direct access to the value of 77, but at the intracellular 
level. Measuring r\ at the scale of several cells either di- 
rectly, or indirectly by measuring E and r and using the 
relation rj — Et, could probably be possible, but we are 
not aware yet of a published method. 

Since £ is a dimensionless geometrical quantity, its 
value is insensitive to cell size: the observed ~ 7- fold in- 
crease in e cannot be related to changes of cell size over 
time. It should involve other causes, that remain to be 
determined. Part of the measured increase in a/rj values 
may be due to cell divisions which decrease the average 
cell area by a factor of ~ 2 between young and old pupae. 
The observed increase is a factor of ~ 25, and thus should 
involve other causes as well. One candidate is a decrease 
of the tissue viscosity r\. However the near constancy of 
the viscoelastic time would then imply a corresponding 
increase in tissue elasticity, a rather unlikely coincidence. 
We therefore expect that stronger tensile forces are at 
work in the tissue in older pupae. From 18 hours APF 
onwards, the lateral part of the scutellum undergoes api- 
cal cell contractions that shape the lateral domain of the 
tissue [29] . This observation led us to divide the ablation 
experiments into the three age groups. It may also ex- 
plain the changes in mechanical properties measured in 
the tissue for middle-aged and old pupae. To precisely 
analyse the causes of the increases in e and a/rj would 
therefore require to identify genes that specifically affect 
the lateral contraction of the tissue. 

Our experimental method should be applicable to 
other tissues. Performing experiments in an epithelium 
which is flat enough to significantly increase L, and with 
a sufficient signal to noise ratio on the estimate of £, 
would enable to check the scaling of £ versus L given by 
Equation 5. 



V. METHODS 
A. Experiments 

Drosophila melanogaster larvae were collected at the 
beginning of metamorphosis (pupa formation) [52]. The 
pupae were kept at 25°C, then dissected and mounted 
as described in [53] (ESM Figure 1). Using the adherens 
junction protein E-Cadherin fused to Green Fluorescent 
Protein [54] we imaged the apical cell junctions by fluo- 
rescent microscopy [53]. 

The time-lapse laser-scanning microscope LSM710 
NLO (Carl Zeiss Microimaging) collected 512x512 pixel 
images in mono-photon mode at 488 nm excitation (pixel 
size in the range 0.24 to 0.29 yum) through a Plan- 
Apochromat 63x1.40 oil objective. Before and after sev- 
ering, the time interval between two consecutive frames 



was either 393 ms or 970 ms according to the scanning 
mode (bi- or mono-directional), without any apparent 
effect on the results presented here. We focused on the 
adherens junctions to be severed; due to the epithelial 
surface curvature especially in old pupae, adherens junc- 
tions of innermost and outermost regions were slightly 
out of focus (Figure 1 and Movies 1-3). 

We defined a region of interest as an annular region 
between two concentric circles. We severed it at the 10 th 
frame (Figure 1) and defined L as the radius of its inner 
circle. Laser severing was performed using Ti: Sapphire 
laser (Mai Tai DeepSee, Spectra Physics) in two-photon 
mode at 890 nm, < 100 fs pulses, 80 MHz repetition rate, 
<~ 0.2 W at the back focal plane, used at slightly less 
than full power to avoid cavitation (ESM Figure 2) [21]. 
Severing itself had a duration ranging from 217 to 1300 
ms (according to the size of the region of interest and the 
scanning mode) during which no image was acquired. We 
could check that the severing had been effective to remove 
their adherens junctions, rather than simply bleach them 
(ESM Figure 2): after retraction the cells at the wound 
margin moved apart more than a cell diameter; moreover, 
there was no fluorescence recovery (Figure 1 (d)). 



B. Analysis 

At the tissue scale, for each frame after the severing the 
inner severed tissue boundary was fitted with an ellipse 
(Figure 1 (c-d)) using the ovuscule [36] Image J plugin 
[55]. We have modified it so that each fitting procedure 
started from the ellipse fitted on the preceding image: 
such tracking improved the speed and robustness. Data 
of Figure 2 are linearly fitted over a 5 frames sliding win- 
dow in order to improve the signal to noise of position 
and velocity determinations. 

Independently, at the cell scale, images were denoised 
with Safir software [56, 57]. We took the logarithm of 
the grey levels, to obtain comparable intensity gradi- 
ents in differently contrasted parts of the image (Fig- 
ure 4 (a)). This allowed us to use Kanade-Lucas-Tomasi 
(KLT) [58, 59] tracking algorithm [60] by selecting - 10 2 
"features" of interest, corresponding to most cell ver- 
tices (Figure 4 (a)). Features were tracked from frame to 
frame: they followed the cell vertex movement towards 
the center, and we have observed no cell neighbour swap- 
ping. We checked that the features' center of mass had 
only a small displacement, which we subtracted from in- 
dividual feature displacements without loss of generality. 
To implement a quasi-lD analysis, we collected the fea- 
tures in a band along the Y axis and similarly along the 
X axis (Figure 4 (b-c)), and both axes were analysed 
separately. The band width (here 1/3 of the initial cir- 
cle diameter) was chosen to be large enough to perform 
large statistics on features, and small enough to treat the 
features displacement as one-dimensional. 
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C. Comparison of model with experimental data i?^ xp (i): 



We numerically solved adimensionalized equations (see 
ESM) using the Matlab solver pdepde. For a given set of 
parameter values (r, £), we calculated Eq from the strain 
of the feature of largest initial amplitude, and rescaled 
the solution of ESM Equation 21 so as to satisfy ESM 
Equation 10. Using the experimental initial positions 
of features R^ Kp (t — + ) and the calculated strain field 
e(z, t), we simulated the trajectories Rf m (t), i = 1 . . . Nf, 
where Nf was the number of features. 



The mismatch £ with our computation was evalu- 
ated by comparison with the experimental trajectories, 



Nf too 
N f l=1 *oc t=1 

where || || denoted the Euclidean norm, and t^, the 
number of movie images. 
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Electronic Supplementary Material (ESM) 

MODEL 

The external dissipation is not accessible to single cell junction ablation experiments and its estimate requires 
modelling [1] . Our model is analogous to models of viscoelastic relaxation proposed for local severing experiments of 
acto- myosin bundles in living cells [2]. 

Hypotheses 

We use an Eulerian description. Since the severed domain area visibly decreases with time, we do not enforce any 
in-plane incompressibility condition. Note also that in-plane incomprcssibility would fix the velocity profile, which 
is incompatible with the experimental observation of different profiles (Figure 5). We rather assume for simplicity 
that the out-of-planc thickness is a passive variable: hence the 3D Poisson ratio is 1/2, but the in-plane effective 2D 
Poisson ratio vanishes. We test this hypothesis by checking a posteriori that the model predictions are compatible with 
experimental observations. Within this simplification, the antero-posterior X and medio-lateral Y axes (Figure 1) 
then decouple. Finally, we define a as the part of stress which has been removed by ablation. 

We thus consider only one spatial dimension, with the variable z denoting position along either X or Y. We model 
the retracting slab of tissue as a linear viscoelastic Kelvin- Voigt material (that is, with a steady response at long time) 
with a Young modulus E and a viscosity r\. We neglect the line tension around the tissue domain: first, because we 
observe that the ellipse does not remain round (this suggests the line tension is weak, without being a proof, since we 
could also imagine that it has the same anisotropy as the strain); second, because even if the line tension exists, it is 
small enough that the ellipse contour remains visibly rough, see Movies 1-3. Two terms contribute explicitly to a: an 
elastic one (which implicitly includes the external pressure) and a viscous one: 

dv 

a(z,t) = Ee(z,t)+r)—(z,t). (8) 
oz 

Taking into account a fluid external friction, with a friction coefficient ( of dimension Pa.s.m -2 , the force balance 
equation reads: 

^(z,t)-Cv(z,t)=0. (9) 
Combining Equations 3, 8 and 9 yields Equation 4. 

Initial, boundary and final conditions 

The initial length of the retracting slab is 2L: z e \—L,L]. Assuming that the initial state respects mechanical 
equilibrium (over the relevant time scale of 10 to 100 s), the initial strain is uniform, and we have for all z: 

e(z,t = 0) = s o , 

a(z,t = 0) = Ee a . (10) 
At the boundary of the severed region, we have for all t: 

a(±L,t) = 0. (11) 
Since the final state is at mechanical equilibrium, we have for all z: 

lim e(z, t) = 0, 

t— >oo 

lim a(z, t) = 0. (12) 

t— >oo 

Injecting Equation 3 into Equation 8, using conditions given by Equations 10-12 yields by integration over time: 

£ (±L,i) = £ e" t/r , (13) 

where the relaxation time is viscoelastic: r = ^. 
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Initial velocity profile 

Equations 3, 10, 13 yield an approximate analytical prediction for the initial velocity profile at time St -C r, as 
follows. An expansion of e at first order in ^: 

e(z,6t)=e +e 1 (z)—, (14) 

T 



defines the initial strain rate as: 

Ei(z) e(z, St) — e de 



(z,t = 0). (15) 



t St dt 

Substituting Equations 14-15 into Equation 4, we obtain an ordinary differential equation for £i(z): 

Dr^-e 1 (z) = 0, (16) 

with a diffusion coefficient D = ^. With the boundary condition e\(z = ±L) = — e , the solution of Equation 16 
reads: 

, , coshfcz 

£lW = - £o ^hfci: (17) 

where k 2 — j^, oi kL = £}/ 2 since £ = Spatial integration of Equation 17 using Equations 3 and 15 yields 

Equation 6. 

Limit regimes 

ESM Figure 4 shows two limit regimes of Equation 4. 

• When internal viscosity dominates: external friction may be neglected, and Equation 9 reduces to ^ = 0, so that 
a is uniform and its value is fixed at the boundaries: <j(z,t) — a(±L,t) — for all (z,t). Integration yields: 

e(z,t)^e a e- t ^. (18) 

The strain is uniform at all times (ESM Figure 4 (A)). Since the center of the slab does not move, v(0,t) — for all 
t, the velocity profile is linear: 

v(z,t) = J' ^(z',t) dz' = z e- t/T . (19) 

As expected, in the limit of small friction £ <C 1, Equation 6 reduces to Equation 19 when t = 0. In this regime, all 
parts of the tissue relax exponentially with the same time scale r. 

• When external friction dominates: the internal viscosity of the tissue may be neglected. Combining Equations 3, 8 
and 9 yields a diffusion equation for the strain field: 

de d 2 e , , 

with conditions given by Equations 10-12. Strain diffuses from the boundaries on a time scale td = L 2 /D = (L 2 /E, 
and the strain field remains inhomogeneous until t>Tc (ESM Figure 4 (B)). 

One may note that the strain diffusion equation (Equation 20) is identical to that obtained in the context of gel 
dynamics [3, 4]. Taking into account interstitial flow (i.e. the transport of water across the epithelium) would only 

modify the diffusion constant: -D pormca tion = E (4 + () , where k is the permeability coefficient. 

Dimensionless variables 

In the general case, to solve Equation 4 numerically we normalize variables by defining: z = z/L, t = t/r and 
£ = e/eq. Equation 4 then becomes: 
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with the initial condition: 

e(z,i=0) = l, (22) 

the boundary conditions: 

e(z = ±l,i) = e-*, (23) 

and the final conditions: 

e(z,f^oo)=0, (24) 

Since £ = td/t, the limit regimes are recovered when t d <C t (internal viscosity dominates, £ <C 1), and t d ^> t 
(external friction dominates, £ » 1). 
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[4] Doi M 2009 Gel dynamics. J Phys Soc Jpn 78:052001. 



14 



ESM FIGURES 




FIG. 1. (ESM) Set-up for live imaging. (A) The pupal case (light brown) is partially removed, while the cuticle (transparent) 
is kept intact. Blue cross: approximate position of severing, along the symmetry axis of the scutellum. (B) The pupa (white) is 
mounted between slide and coverslip with its back facing the objective (dark grey). Antero-posterior axis is horizontal (anterior 
towards the left, posterior towards the right). 





FIG. 2. (ESM) Range of acceptable severing powers. Images shown are just before the severing (top) and 24.3 s after the 
severing (bottom) as a function of the proportion (indicated at top left) of the full laser power, here at 890 nm. Below 20%, 
no severing occurs. From 30 to 50%, severing is partial. From 60 to 80%, severing is complete and cell-cell junctions move 
significantly; experiments presented in this work are performed in this range. Above 90%, cavitation occurs. These ranges vary 
slightly from one experiment to the other, so that the quality of severing must be carefully assessed a posteriori. 
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£, - X axis 



FIG. 3. (ESM) Check of r and £ isotropy. (A) Relaxation time r: for each experiment, represented by one point, we plot the 
value obtained along the Y direction (estimated as in the Inset of Figure 2) versus the value obtained along the X direction. 
Age is color-coded as in Figure 3: green, young; red, middle-aged; blue, old. (B) Same plot for the friction/ viscosity ratio £, 
fitted by Equation 6 as in Figure 5. Note the logarithmic scale. Two values of £ of order 10 -3 and 10 -4 are below the plotted 
range. 
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FIG. 4. (ESM) Limit regimes: space-time profiles of strain fields obtained from numerical simulations of Equation 4 for eo = 1, 
r = 10 s, L = 10 fim. (A) Uniform strain relaxation, when internal viscosity dominates (low f, initial velocity profile similar 
to the blue curve in Figure 5): D = 100 fim 2 .s _l , td = 1 s. (B) Strain diffusion, when external friction dominates (high £, 
initial velocity profile similar to the green curve in Figure 5): D — 1 /im 2 .s" 1 , to = 100 s. Strain value in each grid square is 
color-coded from (dark blue) to 1 (dark red). 
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FIG. 5. (ESM) Error maps in the (r, £) plane, semi-log axes. Experiments are those of Figure 5 and Movies 1-3; error is 
estimated according to Equation 7. The minimal region of the error landscape (its deepest valley) is too wide to allow a direct 
measurement of the parameters. However, depending on the orientation of the minimal region of the error landscape, the 
computed error maps fall broadly into three categories, which agree with the three age groups. There is also a good agreement 
with hatched regions, which are the results of the independent measurements of the parameters r and £ (Figure 3 (c) and ESM 
Figure 3). (A) Small £, old pupa (Movie 3). Minimal valley parallel to the £ axis: this suggests that only an upper bound, of 
the order of 1 (white bar), can be given for the value of £. Friction is negligible. (B) Finite £, middle-aged pupa (Movie 2). 
Minimal valley intermediate, parallel to the r axis: friction cannot be neglected, and £ admits a finite estimate. (C) Large £, 
young pupa (Movie 1). Minimal valley parallel to the r axis, at very large £: this suggests that only a lower bound can be 
given for the values of r and £. 





